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1 Abstract 



A Chow test for structural breaks in the surface temperature series is used 
to investigate two common claims about global warming. Quirk (2009) pro- 
posed that the increase in Australian temperature from 1910 to the present 
was largely confined to a regime-shift in the Pacific Decadal Oscillation 
(PDO) between 1976 and 1979. The test finds a step change in both Aus- 
tralian and global temperature trends in 1978 (HadCRU3GL), and in Aus- 
tralian rainfall in 1982 with flat temperatures before and after. Easter ling 
& Wehner (2009) claimed that singling out the apparent flatness in global 
temperature since 1997 is 'cherry picking' to reinforce an arbitrary point 
of view. On the contrary we find evidence for a significant change in the 
temperature series around 1997, corroborated with evidence of a coincident 
oceanographic regime-shift. We use the trends between these significant 
change points to generate a forecast of future global temperature under spe- 
cific assumptions. 



2 Introduction 



Climatic effects of fluctuations in oceanic regimes are generally studied by 
decomposing rainfall and temperature into periodic components: e.g. sin- 
gular spectrum analysis (SSA) ( |Ghil Vautard 19911, and variations on 
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principle components analysis (PCA) (Parker et al. 2007 1 . Such approaches 
can capture such effects as the influence of short period phenomena like El 
Nino on Australia (Mantua et al. 1997; Guilderson & Schrag 1998), and 



the potential for longer term phenomena such as the Pacific Decadal Oscilla- 



tion (PDO) to 'offset' putative increases in global temperature (Keenlyside 
et al. 2008). Methods designed for finding and testing structural breaks 



address more infrequent regime-shifts. An F-statistic known as the Chow 



test (Chow 1960) based on the reduction in the residual sum of squares 



through adoption of a structural break, relative to an unbroken simple lin- 
ear regression, is a straightforward approach to modeling regime-shifts with 
structural breaks. 



Two claims are evaluated here. In the first, Quirk (2009) proposed a 



model of Australian temperature with a regime-shift (herein Q09): slightly 
increasing to 1976, rapidly increasing to 1979 (the shift), and slowly increas- 
ing ther after. The increase in Australian temperature of around 0.9°C from 



the start of the readily available records in 1910 (BoM 2009) is convention- 



ally modeled as a linear trend and, despite the absence of clear evidence, 
often attributed to increasing concentrations of greenhouse gases (GHGs) 
e.g. (ICSIRO & BoMl 120071). IQuirkl (120091) proposes a different model, with- 



out testing its significance. Regarding only the trends of 0.034°C and 0.05°C 
per decade in the sections before and after the shift substantially lowers the 
underlying rate of warming potentially attributable to anthropogenic global 
warming (AGW). 



Hartmann & Wendler (2005) noted for the climate of Alaska, that recog- 



nition of an abrupt temperature change in climate 1976 profoundly changed 
the underlying temperature trends from strongly positive, to slightly nega- 
tive, contradicting a view of increasing greenhouse gases (GHGs) gradually 
increasing temperatures. The event dubbed the Great Pacific Climate Shift 
(Kerr 1992), is clearly seen in a change in the average polarity of the Pa- 



cific Decadal Oscillation (PDO) from a cool to a warm phase (Hartmann & 



Wendler 2005). 



The second example is a claim of lack of statistical significance. Regard- 
ing the declining temperature since the El Nino event in 1998, Easterling &; 



Wehner (2009) state that using 1998 as a starting point is 'cherry-picking' 



to justify a point of view and therefore, statements such as 'global warm- 
ing stopped in 1997' have no basis. Structural change methods are used to 
determine if the period is a significant change point. 



From a broader perspective, econometrics research into regime-shifting 
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models indicates that detecting breaks in time series is important and ignor- 



ing them can be very costly (Pesaran & Timmermann 2004 ) . The cost in 



climate forecasting would be inaccurate forecasts. That is, while a forecast 
based on a linear model would indicate steadily changing global tempera- 
tures, forecasts based on shifts would reflect the moves to relatively static 
mean values. The choice of underlying model may also impact estimates of 



the magnitude of climate change, as the Quirk (20091 attribution of a trend 



of either 0.04°C or 0.09°C shows. Further, oceanographic changes have been 
postulated as the main causes of changes in global temperature in the last 



30 years (Gray 20071, including the recent flat to slightly declining trend in 



temperatures since 1997. 

A larger relevance of structural change models to global temperature 
data is suggested by the presence of long term persistence (LTP) in com- 
plex natural systems ( Koutsoyiannis & Cohn 2008 Halley| 2009). It is 
known that step changes in the mean at multiple scales can reproduce LTP 
behaviour (Koutsoyiannis 2002 Stockwelll 120061). Change points can be 



found in the mean, or standard deviation (JBaiJ [1994]). 

There are well-known difficulties in forecasting with unstable structural 



breaks models (Hamilton 1989 2007). In particular it is known that a least- 



squares estimator of change points is not reliable for nonstationary 1(d) data 



with 0.5 < d < 1.5 (Hsu & Kuan, 2007). But when there is a change, the 



least squares estimator is known to be reliable on stationary 1(d) data (i.e., 
-0.5 j d j 0.5). Here d is the fractional differencing parameter as introduced 



by Granger & Joyeux (1980). We also have some reservations about the data 



being of adequate length, as models induced from the temperature record 
can be contingent on subsequent data (iStockwell 2009). Given these un- 



certainties, it is important to develop models rigorously, with a transparent 
development sequence satisfying both significance and optimization of the 
model space. 



3 Method 



We attempt to verify the optimality and significance of the breakpoints, 
by reference to both climate data and regime-shifts. Australia- wide annual 
average rainfall and temperature anomalies were downloaded from the BoM 
data site (BoM 2009). Temperature is available as an annual average in 
anomaly °C from 1910 to 2008 while rainfall in anomaly millimeters covers 



1900 to 2008. Quirk (2009) only shows temperature data from 1950 to 2003. 
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We also examine the global monthly temperature anomalies from the Hadley 
Center, HadCRU3GL, for confirmation of the presence of breakpoints at the 



larger scale ( Rayner et al. 2006 1 . Comprehensive evaluation of the full range 



of available global datasets was beyond the scope of this study. 

The approach to developing a structural break model was to calculate the 
F statistic (Chow test statistic) for potential breaks at all potential change 



points using the R package strucchange (Zeileis et al. 2002 1 . In this test 



the error sum of squares (ESS) of a linear regression model is compared 
with the residual sum of squares (RSS) from a model composed of linear 
sections before and after each potential change point. If n is the number of 
observations and k the number of regressors in the model, the F statistic is: 

F = {RSS - ESS)/ESS *{n-2*k) 

Potential breakpoints are indicated by peaks in the value of the F statis- 
tic (dashed vertical lines in figures) . The red line indicates the height above 
which a series with a break at that location would be significantly better (at 
the 95% confidence limit) than a single linear regression. 



Table 1: Fit of models as described in text and shown in Figures la,b and 2a,b. R2 
is the regression correlation coefficient, R2b is the correlation coefficient with one 
break, Pb is significance of improvement of the break, over the linear model, PI and 
P2 and the significances of the slopes of the first and second segments, respectively. 





R2 


Breakdate 


R2b 


Pb 


PI 


P2 


Aus. Temperature 


0.384 


1978.00 


0.482 


0.000 


0.151 


0.139 


Aus. Rainfall 


0.041 


1982.00 


0.164 


0.017 


0.877 


0.026 


Global Temp. 1910-2008 


0.721 


1978.00 


0.785 


0.000 


0.000 


0.000 


Global Temp. 1976-2008 


0.613 


1997.08 


0.672 


0.000 


0.000 


0.409 



The relative fit of the models is indicated with a regression coefficient 
i? 2 , as listed in Table [1} In addition, p values indicate the significance of 
model improvements (Pb), and the non-zero slope parameters (PI and P2). 

While the strucchange software also implements a search for multiple 
breakpoints, we adopt a conservative approach without undue complexity, 
mindful of the problematic nature of simultaneous validation of multiple 
structural break models (Zeileis &; Kleiber, 2004 1 . The optimal position of 
single breaks is clearly indicated with a plot of the F statistic aligned with 
the x-axis of a plot of the time series. 
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4 Results 



As a preliminary check, the estimates of the fractional differencing parameter 
d on the Australian temperature and rainfall data were estimated at 0.31 
and 0.13 respectively (using R package fracdif /). Therefore, the Chow 
test is expected to be reliable. The estimate of d on HadCRU3GL global 
temperature series is 0.50, so results must be interpreted more cautiously. 

The F statistics, breakpoints, confidence intervals, the simple regression 
and segmented models for annual Australian temperature from 1910 to 2008, 
are shown on Figure [T^,. The peak in the F statistic above the red line indi- 
cates that 1978 is the optimal location for a break, indicated with a dashed 
vertical line. The 95% confidence interval spans 1974 to 1981, indicated with 
a red bracket on the x-axis of the F statistic. The structural break model 
is a significant (p < 0.001) improvement in the fit of the model as shown 
by the increase in R 2 from 0.38 to 0.48. The slope of the two segments is 
not significant (p = 0.15 and p = 0.14). Because of the significance of the 
break, and the non-significence of slopes in the segments, a model with a 
single step, or structural change model is justified (Figure [T^ blue). 

Regarding the precipitation series, a single break-date at 1982 (Fig- 
ure lip) significantly improves the fit (p = 0.02) as indicated by the increase 
in Ft from 0.041 to 0.164 (Table [TJ. The trend in the segment before the 
break is not significant {p = 0.88), but a decline in precipitation after the 
break is significant {p = 0.03), most probably due to the anomalously high 
rainfall in the years after 1978. 

A broad peak in the F statistic for annual global surface temperature 
from 1910 at 1978 (Figure [2^,) (p < 0.001) improves the fit of the model 
from an R 2 of 0.56 to 0.68. The 95% confidence interval spans 1966 to 
1985, indicated with a red bracket below. When the monthly global surface 
temperature series is truncated to the break-date of 1978, the F statistics 
indicate a strong break-date in 1997(8) (Figure [2p), R 2 increasing from 
0.61 to 0.67) with a rising trend prior to 1997, and a fiat trend thereafter 
(p = 0.41). The 95% CI spans 1996(8) to 1997(12). 



5 Discussion 

These results confirm a break-date around 1978, as suggested by the general 



form of Q09 (Quirk 2009 1, but differs in respect to one, not two breaks, with 
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1910 1920 1930 1940 1950 1960 1970 1980 1990 2000 2010 



Figure 1: Breakpoints in Australian temperature (a) and rainfall (b). The upper 
section of each figure is the F statistic from the Chow test for breaks with optimal 
breaks (vertical, dashed), level of significance (red line) and confidence interval (red 
bracket). The Australian temperature (below), with the OLS regression model 
(green) and the optimal break model (blue). 
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non-significant trends outside of the optimal break-date. This is known as 
a change point model, characterized by abrupt changes in the mean value 
of the underlying dynamical system, rather than a smoothly increasing or 
decreasing trend. The confidence in 1978 as a break-date is further strength- 
ened by the results for global temperatures since 1910, indicating the se- 
ries could be described as gradually increasing to 1978 (0.05 ± 0.015°C per 
decade), with a steeper trend thereafter (0.15 ± 0.04°C per decade). 

The Chow test since 1978 finds another significant break-date in 1997, 
delineating an increasing trend up to 1997 (0.13 ± 0.02°C per decade) and 
non-significant trend thereafter (—0.02 =L 0.05°C per decade). Contrary to 
claims in Easterling & Wehner ( |2009 ) that the 10 year trend since 1998 is 
arbitrary, structural change methods indicate that 1997 was a statistically 
defensible beginning of a new, and apparently stable regime. 

The significance of the dates around 1978 and 1997 to climatic regime- 
shifts is not in dispute, as they are associated with a range of oceanic, 
atmospheric and climatic events, whereby thermocline depth anomalies as- 
sociated with PDO phase shift and ENSO were transmitted globally via 



ocean currents, winds, Rossby and Kelvin waves (Guilderson & Schrag 1998 



McPhaden k Zhang 2004 Wainwright et al. 2008). Even though the pat- 



tern of response would vary in different parts of the globe, a step-change in 



temperature remarkably similar to Australia also occurred in Alaska (Hart 



mann & Wendler, 2005 1 . 



Given the assumption that AGW is not a cause of regime-shifts, an 
assumption that would require further verification, these results suggest a 
considerably stronger claim for Australian climate, i.e. because the slope of 
the two segments is not significant (p=0.15 and p=0.14), any trend due to 
increasing CO2 is statistically insignificant. 

Decomposition of global temperature into regime-shifts and AGW is 
more ambiguous, but it would be argued, similar to Q09, that a large propor- 
tion of the increase in temperature increase between 1978 and 1997 was asso- 
ciated with the propagation of a major regime-shift, with low rates of AGW 
in the segments before and after. Supporting this view are such oceano- 
graphic observations as the large reduction in the rate of Pacific equatorial 
upwelling around 1976 and its resumption in 1998, a reduction of westward 
volume transport between Australia and Indonesia of 23% from 1976-7 and a 



general elevation of sea temperature by 1-2°C (Wainwright et al. , 2008). Fol- 



lowing two decades of weaker flow coinciding with rising global temperatures 
from the break-date of 1978, convergence of cold interior ocean pycnocline 
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water towards the equator increased from 13.4 to 24.1m 3 s , and Pacific 



sea surface temperatures cooled (McPhaden & Zhang, 2004). Cai et al 



(2007) asserts a role for strengthening of the global conveyor in global tem- 



peratures, by intensifying the rate of heat transfer out of the off-equatorial 
region and into the subtropics. The relative contribution and interaction 
of various ocean phases: IPO, PDO and ENSO, the duration of lags, and 



precipitation effects are more uncertain (Power et al. 1999). 



One can still argue for human influence on regime-shifts. Cai et al. 



(2007) attributes changes in the global conveyor to anthropogenic aerosols, 



and Vecchi et al. (2006) attributes a weakening of the Walker circulation to 



anthropogenic forcing. Neither attribute abrupt regional climate changes to 
non-natural causes. There is also the view that stable temperatures since 
1997 may be a result of an 'offset' or masking of AGW effect by natural vari- 



ation (Keenlyside et al. 2008 ). For consistency, this view must also entertain 



the possibility of putative AGW being amplified by natural variability in the 
previous regime from 1978 to 1997. 

Assuming a regime-shift from 1978-98 reduces the estimates of the under- 
lying rate of AGW warming from around 0.14°C to 0.05°C per decade. An 
increase of 0.5°C by 2100 is consistent with low-end empirical estimates of 



climate sensitivity, such as Spencer & Braswell (2008) at 0.6°C for IXCO2 



but considerably lower than the IPCC projections for the most common 
AGW scenarios diPCCl[2007]. Our model is consistent in period and timing 



to the two component model proposed by Akasofu (2009), with a contribu- 



tion from a natural multidecadal oscillation of 0.15°C per decade between 
1975 and 2000, and a steady contribution of 0.05°C per decade from the 
upswing of a multi-centennial oscillation, such as the natural cycle of glacial 
formation and regeneration called the Great Season Climatic Oscillation 



(Boucenna 2008). In our case, the model was developed from statistically 



justified, empirical analysis of temperature data. 

Figure [3] illustrates the prediction for temperatures to 2100 following 
from our structural break model, the assumptions of continuous underlying 
warming, regime-shift from 1978 to 1997, and no additional major regime- 
shift. The projections formed by the trend to 1978 (presumed AGW warm- 
ing, green) and the trend in the current regime (AGW offset by regime-shift, 
cyan) predicts constant temperatures for fifty years to around 2050, similar 
to the period of fiat temperatures from 1930-80, then increasing to about 
0.2°C above present by 2100. It must be kept in mind that this extrapola- 
tion is based on greatly simplified assumptions regarding regime-shifts and 
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1950 2000 2050 2100 

Figure 3: Prediction of global temperature to 2100, by projecting the trends of 
segments delineated by significant regime-shifts. The flat trend in the tempera- 
ture of the current climate-regime (cyan) breaks upwards around 2050 on meeting 
the (presumed) underlying AGW warming (green), and increases slightly to about 
0.2°C above present levels by 2100. The 95% CI for the trend uncertainty is dashed. 
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trends, and does not incorporate many of the complexities and natural forc- 
ing factors operating to produce climate change in the real world (Marsh 
20091). 



These results justify further development of more complex regime-shift 
models of temperature and rainfall for forecasting purposes, including at- 
tempting to decompose global climate into temporally and spatially differ- 
entiated regime-shift models. 
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